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MHD WAVES IN CORONAL HOLES 


1 2 
D. Banerjee , S. Krishna Prasad 


Coronal holes are the dark patches in the solar 
corona associated with relatively cool, less dense 
plasma and unipolar fields. The fast component of 
the solar wind emanates from these regions. Several 
observations reveal the presence of magnetohydrody¬ 
namic (MHD) waves in coronal holes which are be¬ 
lieved to play a key role in the acceleration of fast 
solar wind. The recent advent of high-resolution 
instruments had brought us many new insights on 


the properties of MHD waves in coronal holes which 
are reviewed in this article. The advances made in 
the identification of compressive slow MHD waves in 
both polar and equatorial coronal holes, their pos¬ 
sible connection with the recently discovered high¬ 
speed quasi-periodic upflows, their dissipation, and 
the detection of damping in Alfven waves from the 
spectral line width variation are discussed in partic¬ 
ular. 


1. INTRODUCTION 

Solar corona, when seen in ultraviolet (UV) and 
X-ray radiation, show dark patches on the Sun, 
called coronal holes. They appear dark only in 
the layers of the solar atmosphere with tempera¬ 
tures exceeding 10^ K, making them indistinguish¬ 
able from surroundings in the photospheric and low 
chromospheric images and hence the name ‘coro¬ 
nal’ holes. These regions are cooler and less dense 
than their surroundings and are associated with 
open magnetic fields of the Sun. Magnetic flux is 
mainly unipolar with huge imbalance between the 
positive and negative fluxes [Wang^ 2009]. This is 
connected with the migration of decaying active re¬ 
gion fields towards poles, and consequently, the lo¬ 
cation of coronal holes depends on solar cycle. Dur¬ 
ing solar minimum, coronal holes are mostly con¬ 
fined to polar regions (polar coronal holes) while 
they can be found at lower latitudes (equatorial 
coronal holes) as well, during solar maximum. The 
fast component of the solar wind emanates from 
the coronal hole regions. Hence, their presence in¬ 
fluences the speed of slow solar wind coming from 
equatorial regions which reaches up to 600 km s“^ 
[Zhang et al.^ 2003] during solar maximum, faster 
than the average speeds ( 400 km s ^) [Woch 
et al.^ 1997]. Interested readers are referred to an 
excellent review by Cranmer [2009], for further in¬ 
formation on coronal holes. 
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Polar coronal holes are large and persistent for 
most part of the solar cycle. Solar wind emanating 
from these regions is fast (^^750 km s“^), tenu¬ 
ous and relatively homogeneous compared to the 
slower, denser and highly variable wind emanating 
from low latitude regions [Woch et al.^ 1997; Mc- 
Comas et a/., 2000]. During solar minimum, the 
global structure of solar wind is more or less uni¬ 
form with the highly variable slow wind restricted 
to a narrow belt around equatorial regions whereas 
during solar maximum, it becomes complex with 
flows arising from multiple sources (streamers, ac¬ 
tive regions etc.,) at all latitudes [McComas et al.^ 
2002, 2008]. 

One of the striking features visible in the polar 
coronal holes is bright ray like structures above the 
limb called polar plumes that were first noticed in 
white light eclipse observations [van de Hulst^ 1950; 
Saito^ 1958, 1965]. With predominantly unipolar 
magnetic flux concentrations at their foot points 
[Deforest et al.^ 1997], plumes are believed to trace 
open magnetic helds in solar corona [Banaszkiewicz 
et a/., 1998; Young et a/., 1999]. The dark re¬ 
gions between plumes are called interplume re¬ 
gions. Plumes are cooler but denser than the sur¬ 
rounding interplume regions [Wilhelm^ 2006]. It 
is often debated whether plumes or interplumes 
are preferred sites for the acceleration of fast so¬ 
lar wind. While some researchers [Walker et a/., 
1988; Gabriel et a/., 2003, 2005] believe plumes to 
be a major source of high speed wind, magnetic 
flux observations at the foot points [Wang et al.^ 
1997] combined with lower outflow speeds [ Wilhelm 
et al.^ 2000; Patsourakos and Vial^ 2000] and nar¬ 
rower spectral prohles In plumes [Antonueei et al.^ 
1997; Banerjee et a/., 2000a; Teriaea et a/., 2003] 
tend to support otherwise. The reason for this 
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ambiguity can be partly attributed to our lack of 
knowledge on the exact physical process behind the 
acceleration of fast solar wind. Magnetohydrody¬ 
namic (MHD) waves are believed to be one of the 
possible candidates for transporting and deposit¬ 
ing the required energy for acceleration. In coro¬ 
nal holes, plumes are regions of Alfven speed min¬ 
ima (a consequence of the pressure balance across 
these dense structures) which makes them natu¬ 
ral guides for different MHD waves [Nakariakov^ 
2006]. Observations of these waves had been revo¬ 
lutionized by the launch of Solar and Heliospheric 
Observatory [SOHO, Domingo et a/., 1995] in 1995. 
Since then, MHD waves were detected by several 
authors both from remote sensing and in situ ob¬ 
servations. Banerjee et al [2007, 2011] provide a 
comprehensive review of these observations. The¬ 
oretical models were also developed to understand 
the actual energy transfer process by these waves 
[see the review by O/man, 2005]. 

Observations of wave activity, relatively sim¬ 
ple geometry, and their connection to the fast 
solar wind, makes coronal holes interesting to 
study. In this article, we review the observations 
of MHD waves in coronal holes focusing primar¬ 
ily on progress made in the past few years. We 
also discuss on the new evidences of wave damping 
and highlight how such observations can be used to 
probe coronal conditions through seismology. Seis¬ 
mology is an interesting application of observing 
waves in the solar atmosphere but the subject is 
beyond the scope of the present review. 

2. WAVES IN POLAR CORONAL HOLES 

During solar minimum, coronal holes are mostly 
confined to the polar regions. The off-limb part of 
the polar coronal holes appear structured with the 
presence of plume and interplume regions while the 
on-disk part is mostly dark with a few small bright 
patches (coronal bright points) appearing here and 
there. Numerous observations using imaging and 
spectroscopic techniques, revealed the presence of 
different MHD waves in these structures that can 
be categorized into compressive and incompressive 
waves. 

2.1. Compressive Waves 

Compressive waves, as the name suggests, are 
perturbations that cause ffuctuations in density, 
which thereby alter the observed emission or in¬ 
tensity. First such observations in polar plumes 
were made by Withbroe [1983] who found short 
period (a few minutes) variations in Mg X emis¬ 
sion observed by the Harvard Skylab experiment. 
However, the author did not find any correspond¬ 
ing variations in O VI line which is equally den¬ 


sity sensitive and it was suggested that if these 
variations in Mg X emission are real they could 
be due to temperature ffuctuations in response to 
the variations in local plasma heating rate. Af¬ 
ter the launch of SOHO, Ofman et al [1997] ob¬ 
served ffuctuations in polarized brightness (a mea¬ 
sure of electron density) over a coronal hole high 
above the limb between 1.9 and 2.45 i?Q, using 
data obtained with the white light channel of Ul¬ 
traviolet Coronagraph Spectrometer (UVCS). This 
was the first undoubted detection of compressive 
waves in polar coronal holes. Deforest and Gur- 
man [1998] studied the time evolution of several 
plumes by constructing evolution charts from a se¬ 
quence of images observed by Extreme ultravio¬ 
let Imaging Telescope [EIT, Delaboudiniere et a/., 
1995] onboard SOHO. Evolution charts were made 
by taking a strip along a plume at each instant 
and stacking them in time. They found diagonal 
stripes in these charts representing the propagation 
of enhancements and suppressions in brightness 
and suggested them as compressive waves propa¬ 
gating along the plumes. The period at which the 


AIA 171 A 20-July-2010 02:00 UT 



Time (min) 


Figure 1. Top: Locations of all the slits chosen over 
several plume/interplume regions seen at a south polar 
region. Contours of different intensity levels overplot- 
ted as dotted curves, show the exact locations of plumes 
and interplumes. Slits over plume regions are 36'' wide 
and those over interplume regions are 18" wide. Bot¬ 
tom: Time-distance maps constructed from slits 2 & 9 as 
annotated. Adapted from Krishna Prasad et al. [2011]. 
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stripes are repeated is 10 - 15 min and the propa¬ 
gation speeds obtained from the inclination of the 
stripes, were 75 - 150 km s“^. The amplitudes 
were found to be 10 - 20% of the overall intensity. 
Further analysis indicates an increase in wave am¬ 
plitude with height. In a gravitationally stratified 
medium, the falling density will naturally cause an 
increase in wave amplitude if the energy flux car¬ 
ried by the waves is conserved. The observed in¬ 
crease is in agreement with the MHD simulations 
for propagating slow waves leading to the inter¬ 
pretation of these waves as slow magneto-acoustic 
waves [Ofman et a/., 1999]. A number of studies 
followed reporting such oscillations in plumes, in¬ 
terplumes, and coronal holes, using spectroscopic 
data obtained with Coronal Diagnostic Spectrom¬ 
eter (CDS) and Solar Ultraviolet Measurements of 
Emitted Radiation (SUMER) onboard SOHO and 
Extreme ultraviolet Imaging Spectrometer (EIS) 
onboard Hinode [Banerjee et a/., 2000b, 2001a, b; 
OBhea et a/., 2006, 2007; Banerjee et a/., 2009a; 
Gupta et al.^ 2009]. 

Using Extreme UltraViolet (EUV) observations 
from Solar TErrestrial RElations Observatory 
(STEREO), Meintosh et al [2010] analysed sev¬ 
eral plume structures and suggested that prop¬ 
agating disturbances (PDs) along polar plumes 
could be due to the collimated high-speed quasi- 
periodic plasma jets that have similar properties 
rather than due to slow magneto-acoustic waves. 
Further, they conjectured that these jets could 
be responsible for loading a significant amount of 
heated plasma into the fast solar wind. It is im¬ 
portant to clearly understand the nature of these 
PDs to find their contribution (mass/energy load¬ 
ing) to the solar wind. In an attempt to re¬ 
solve this Krishna Prasad et al [2011] studied 
the properties of PDs along several plume and 
interplume regions using data obtained with At¬ 
mospheric Imaging Assembly (AIA) onboard So¬ 
lar Dynamics Observatory (SDO). The authors 
followed the regular time-distance analysis, how¬ 
ever, they chose wider (several tens of AIA pixels) 
slits along plume/interplume regions and averaged 
across them while constructing the time-distance 
maps. Figure 1 displays a south polar region with 
chosen slits over several plume and interplume re¬ 
gions marked (top panel). The width of the slits 
shown are 60 AIA pixels (~ 36^^) for plume region 
and 30 AIA pixels (~ 18'') for interplume region. 
The bottom two panels of the figure show time- 
distance maps corresponding to slit 2 (plume) and 
slit 9 (interplume). PDs can be clearly seen in these 
maps and were found to be unaffected by variations 
in the slit width. The faint plasma jets which oc¬ 
cur at random should not be present in these maps 
as the spatial average across the slit makes them 
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Figure 2. Enhanced time-distance maps for all the pa¬ 
rameters of Ne VIII line and for line width of N IV line as 
annotated. The slanted solid lines marked along the PDs 
indicate a periodicity of 14.5 min and propagation speed 
of 60 km s“^. Ellipses in the bottom panel show periodic 
enhancements in line width at the base of the PDs. Erom 
Gupta et al. [2012]. 

fainter. On the other hand, the stronger jets were 
identified visually and eliminated. The fact that 
the observed PDs are unaffected by variations in 
slit width implies a coherent behaviour. The pro¬ 
jected propagation speeds were in the range of 100 
- 170 km s“^ and were found to be temperature 
dependent. Also, the speeds were slightly higher 
in the interplume regions possibly due to relatively 
hotter temperatures. Based on these properties, 
the authors infer that the observed PDs are due 
to slow MHD waves. Further, this study shows a 
way to avoid the contamination from plasma jets 
by using wider slits. 

Gupta et al [2012] observed PDs in an on-disk 
part of a south polar coronal hole using spectro¬ 
scopic data obtained with SUMER. The spectral 
lines Ne Vlll 770 A, and N IV 765 A, which form at 
lower corona and transition region heights respec¬ 
tively, were used in this analysis. Time-distance 
maps for Ne Vlll show clear signatures of PDs in 
intensity and Doppler shift but are barely visible 
in line width. Similar maps for N IV line did not 
show any propagating feature in any of the line pa¬ 
rameters. However, the map corresponding to line 
width shows periodic enhancements at the bottom 
of the PDs observed in Ne Vlll line. Figure 2 dis¬ 
plays the time-distance maps of all the three line 
parameters for Ne Vlll line and that of line width 
for N IV line. The slanted solid lines indicate a 
propagation speed of 60 ± 4.8 km s“^ and a period¬ 
icity of 14.5 min. The ellipses in the bottom panel 
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mark the periodic increase in line width. Note the 
foot point is at the top with limb towards the bot¬ 
tom so the slanted ridges indicate PDs propagating 
outward. This was the first time PDs were simulta¬ 
neously observed in intensity and Doppler shift in a 
polar coronal hole. The enhancements in intensity 
correspond to a blue shift in velocity which is ex¬ 
pected for an outward propagating slow magneto¬ 
acoustic wave. Further, line profiles were found 
to be pretty symmetric with very small asymme¬ 
tries that were shown to arise from the finite sam¬ 
pling and noise in the data. The authors conclude 
that the line width enhancements at the bottom 
of the PDs could be possibly due to some unre¬ 
solved explosive events while the PDs at greater 
heights show dominant wave signature similar to 
the case observed by Nishizuka and Hara [2011] 
for an active region. A 3D MHD model developed 
by Ofman et al [2012] for a bipolar active region, 
also supports these results. 

Wang et al [2012] studied PDs in coronal loops 
using Hinode/EIS observations and forward mod¬ 
els. They analysed the excess emission in the line 
profiles by subtracting a background profile instead 
of the regularly used Red-Blue (RB) asymmetry 
analysis method. Their results favor the wave in¬ 
terpretation for the observed PDs. From a statis¬ 
tical study, Tian et al [2012] found two types of 
oscillations, one close to the loop foot points which 
show coherent oscillations in all the line parame¬ 
ters with blueward asymmetries in the line pro¬ 
files, and the other associated with upper parts 
of the loops which show prominent oscillations in 
Doppler shift, but no significant variation in in¬ 
tensity and line width and no obvious asymme¬ 
tries in the line profiles. The authors interpret the 
former type as due to episodic high-speed upflows 
and the latter type as kink/Alfven waves. Wang 
et al [2013] performed 3-D MHD simulations of 
PDs in long fan loops, and modeled a section of 
a loop with open boundary conditions at the top. 
They injected repetitive tiny flow pulses at the base 
of the coronal loop, with an energy frequency dis¬ 
tribution that follows the flare power-law scaling. 
It was found that the injected upflow pulses in¬ 
evitably excite slow magneto-acoustic waves in the 
loop. The upflows were found to rapidly decelerate 
with height and the resulting PDs are mostly dom¬ 
inated by the wave signatures except near the base 
where the contribution from flows can be notable. 
Based on these results, the authors suggested that 
nanoflare-like impulsive heating events at the base 
of the loop may produce both upflows and waves. 

Following a different approach, Su et al [2014] 
studied the nature of PDs observed at both north 
and south polar regions using 21 datasets. By fit¬ 
ting an exponentially growing sine curve to the 
spatial variation of PDs very close to the limb (0- 


9 Mm), they measured the pressure scale height 
and wavelength in three AIA channels 171 A, 
193 A, and 211 A. The temporal variation was 
then used to estimate the oscillation period in these 
channels. Average values of scale height obtained 
for all the sets of data are 6.4±2.8, 4.8±3.2, and 
3.3±1.7 Mm, respectively in 171 A, 193 A, and 
211 A channels. Corresponding wavelength val¬ 
ues are 35±7, 26±11, 25±20 Mm. The periodicity 
on the other hand, is roughly same 24 min) in 
all three channels. Although the uncertainties in 
the scale height and wavelength are large, the au¬ 
thors had statistically shown that they are actually 
temperature dependent wheras the periodicity re¬ 
mains independent of temperature. Further, they 
calculated the propagation speeds from the ratios 
of wavelength to time period and compared those 
with the independent estimates obtained from the 
pressure scale height (through its relation to tem¬ 
perature). A good agreement was found between 
these values. Depending on these properties, the 
authors suggest that the detected PDs are slow 
magneto-acoustic waves in nature. However, the 
opposite (to expected) dependence of these param¬ 
eters on temperature, i.e., shorter scale heights and 
wavelengths in hotter channels, is not explained by 
the authors. 

[Gupta et al.^ 2010] found accelerating PDs in 
interplume regions. The authors combined spec¬ 
troscopic observations from SUMER and 40" slot 
images from EIS to study the properties of PDs 
in plume and interplume regions. Their analy¬ 
sis reveals PDs of 15-20 min periodicity in inter¬ 
plume regions travelling with projected speeds in¬ 
creasing from 25±1.3 km at their origin to 
330±140 km s“^ at around 160" above the limb. 
These speeds were also found to be temperature 
independent. PDs of periodicity in the same range 
were found in the adjacent plume region but with 
no significant acceleration. Although the observa¬ 
tions were mainly in intensity, the acceleration and 
other properties (such as the supersonic speeds and 
their temperature independence) led the authors 
to conclude the PDs observed in interplume region 
as due to either Alfvenic or fast magneto-acoustic 
waves and those in the plume region as due to slow 
magneto-acoustic waves. 

Slow magneto-acoustic waves were also observed 
in the on-disk part of a coronal hole [O^Shea et al.^ 
2007; Gupta et a/., 2009]. Spectroscopic observa¬ 
tions in N IV 765 A(transition region) and Ne Vlll 
770 A(low corona) lines, obtained with SUMER, 
were used by Gupta et al [2009] to study the phase 
delays between intensity and velocity oscillations. 
They found upward and downward propagating 
waves in the internetwork regions and only upward 
propagating waves in the network regions of the 
polar coronal hole. 
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2.2. Incompressive Waves 

Incompressive waves do not cause any pertur¬ 
bations in intensity, instead they alter other spec¬ 
tral parameters. Therefore, spectroscopic observa¬ 
tions are necessary to observe them. For instance, 
the transverse Alfven waves which are incompress¬ 
ible in the linear limit, cause oscillations only in 
Doppler shift. However, the emission in polar coro¬ 
nal holes, being low, often require longer exposures, 
spatial, and temporal averaging to obtain a good 
signal. This makes the Alfven waves unresolved (if 
the time period is shorter than the temporal bin) 
and increases Doppler width rather than causing 
a Doppler shift. In the WKB approximation, the 



Distance from Sun Center (Rq) 


< 

I 


jj) 

o 

D) 

k. 

3. 

X 

_D 

LL 

C 

CD 

< 





U . I . I ......... 1 ..... . . 7 ■•■a-L...... ■ ■ ■ ■ I 

0.0 0.1 0.2 0.3 0.4 0.5 


h(R) 

o 


Figure 3. Top: Variation of line width (left) and non- 
thermal velocity (right) with height above the limb for 
the spectral line Ee Xll 195.12 A. Solid and dotted lines 
represent values before and after the subtraction of Ee Xll 
stray light profile, respectively. Dashed lines represent 
parabolic fits to the corrected data. Note the non-thermal 
velocities shown here are devoid of instrumental broad¬ 
ening but the line width values are not. Bottom: Alfven 
wave energy flux as a function of height above the limb. 
Values shown in different line styles correspond to dif¬ 
ferent non-thermal velocities shown in the top panel. 
The calculations were done for a constant magnetic field, 
G. Adapted from Bemporad and Abbo [2012]. 


Figure 4. Variation of Doppler width (after instrumen¬ 
tal correction) as a function of radial distance for sev¬ 
eral spectral lines as indicated in the plot. The dashed 
lines represent the corresponding values expected for un¬ 
damped wave propagation. From Hahn et al. [2012]. 


Alfven wave energy flux [F) crossing a surface of 
area A, can be written as 

F = p{Sv^)vaA = .[^{5v^)BA (1) 

V Att 

where va = B/^Airp is Alfven wave velocity, p is 
plasma mass density, B is field strength, and {5v‘^) 
is the mean square velocity related to the observed 
non-thermal velocity as = {6v‘^)/2. It may be 
noted that this approximation is not well applica¬ 
ble when wave reflections or refractions are strong, 
and when damping is important. Nevertheless, it 
is useful to describe the general propagation. For 
a flux tube geometry BA remains constant, and 
if the energy flux is conserved (i.e., if the Alfven 
waves are undamped as they propagate), then the 
non-thermal velocity ^ follows 

(2) 

Equation 2 implies, the observed Doppler width 
due to (unresolved) Alfven waves increases with 
fall in density as they propagate outwards. This 
particular feature turned out to be an important 
(although indirect) evidence for the presence of 
Alfven waves. Several authors observed an increase 
of line width with height in polar coronal holes 
suggesting the interpretation of outward propagat¬ 
ing undamped Alfven waves [Hassler and Moran^ 
1994; Banerjee et a/., 1998; Wilhelm et a/., 2004; 
Banerjee et a/., 2009b]. A few others observed an 
initial increase followed by a decrease or level off at 
greater heights [O’Shea et a/., 2003, 2005; Mierla 
et al.^ 2008] and indicated that the decrease in line 
width could be due to Alfven wave damping. 

O’Shea et al. [2005] also observed a turnover in 
the line width variation and found a change in ex- 
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citation mechanism from collisional to radiative at 
the turning point. Singh et al [2003] found in¬ 
crease and decrease of line widths in different lines 
of a simultaneously observed pair and suggested 
that the variation in line width may always not be 
related to propagating waves. Dolla and Solomon 
[2008] found that the decrease in line width is pri¬ 
marily due to scattered light that was not prop¬ 
erly estimated. Further, the increase in line widths 
can also be caused by increase in ion temperatures 
[Seely et al.^ 1997; Tu et a/., 1998]. All these stud¬ 
ies made the detection of Alfven waves in coronal 
holes inconclusive. 

Recently, Bemporad and Abbo [ 2012 ] studied the 
variation of line width in a polar coronal hole us¬ 
ing spectroscopic data from EIS. The observations 
were made with 1 " slit in sit-and-stare mode for 
nearly 21.6 hrs. The exposure times were kept long 
(500 s) to gather the signal up to farther off-limb 
distances. After careful subtraction of the stray 
light estimated using a transition region line, the 
authors found that the width of the Fe Xll 195.12 A 
line increases up to 2^0.14 Rq and decreases higher 
up. Assuming ionization equilibrium, they also cal¬ 
culated non-thermal velocity as a function of height 
which shows similar behaviour. Figure 3 displays 
the variation of line width and non-thermal veloc¬ 
ity with height above the limb for Fe xll 195.12 A 
line (top panel). The solid and dotted lines in the 
plots correspond to the results before and after the 
subtraction of additional stray light caused by a 
blend in the Fe xll line, respectively. Dashed lines 
represent parabolic fits to the corrected data. The 
line width values shown here, were not corrected 
for the instrumental broadening, although the cor¬ 
rection was done before the computation of non- 
thermal velocities. These plots clearly show the 
decrease of non-thermal velocity after an initial in¬ 
crease. It should be noted that the non-thermal 
velocity is computed here by assuming a constant 
ion temperature at all altitudes, so any possible 
increase in ion temperature with height will fur¬ 
ther reduce the non-thermal velocities making the 
decrease even steeper. The authors interpreted 
this behaviour as due to damping in Alfven waves 
and estimated the energy flux assuming a constant 
magnetic field B=8 G. The computed values as a 
function of height were shown in the bottom panel 
of Figure 3. Values shown in different line styles 
correspond to different non-thermal velocities plot¬ 
ted in the top panel. Apparently, the energy flux 
continuously decreases even in the low altitude re¬ 
gion where the line width increases with height. 
It goes down from about ?^1.2xl0^ erg cm“^ s“^ 
at 0.03 i ?0 above the limb to about ^^8.5x10^ erg 
cm“^ s“^ at 0.4 Rq above the limb. In an in¬ 
dependent study, Hahn et al [2012] also found a 
decrease in line width with height following an ini¬ 


tial increase in a polar coronal hole. The scattered 
light contribution was carefully removed and the 
authors rule out any effects due to line of sight 
and changes in excitation mechanism suggesting 
the Alfven wave damping. Figure 4 displays their 
results for several spectral lines. The turning point 
seems to be at different heights for different lines 
indicating multiple structures in the line of sight. 
Dashed lines in the figure mark the respective val¬ 
ues for undamped wave propagation. Assuming 
thermodynamic equilibrium, the energy flux esti¬ 
mated was about ?^7.0xl0^ erg cm“^ s“^ at 0.1 Rq 
above the limb to about ^^ 1 . 1 x 10 ^ erg cm“^ s“^ 
at 0.3 i ?0 above the limb. The dissipated energy 
between these heights amounts to 70% of the en¬ 
ergy required to heat the coronal hole and accel¬ 
erate the fast solar wind (^^ 8 x 10 ^ erg cm“^ s“^; 
Withbroe and Noyes [1977]). However, the au¬ 
thors cautioned that the estimated values would 
be upper bounds since the actual ion temperatures 
can be higher than the electron temperatures used 
here for the estimation of non-thermal velocity. Al¬ 
though suffered by some observational uncertain¬ 
ties and involve certain assumptions, both these 
studies seem to provide a compelling evidence for 
the existence and damping of Alfven waves in coro¬ 
nal holes. 

3. WAVES IN EQUATORIAL CORONAL 
HOLES 

Coronal holes are also observed in the equato¬ 
rial regions during the maximum phase of the so¬ 
lar cycle. Equatorial coronal holes are usually as¬ 
sociated with the remnant active regions. Their 
smaller spatial extent makes them easy to observe 
when they are at the disk center, which however, 
poses a challenge for observations of propagating 
waves expected for an open field configuration. In- 
compressive waves with perturbations perpendic¬ 
ular to the line of sight are difficult to observe, 
nevertheless, there are a few studies which report 
the presence of compressive waves. O’Shea et al. 
[2007] performed a statistical study to find evi¬ 
dence for waves in equatorial coronal holes using 
observations from CDS/SOHO. Oscillations in in¬ 
tensity and Doppler shift were studied for the spec¬ 
tral lines O V 629 A, Mg x 624 A, and Si xii 520 A. 
By measuring phase delays between different line 
pairs, propagation speeds were estimated to be in 
the range of 50 to 120 km s“^. It was also observed 
that the oscillations occur preferentially in bright 
regions that are associated with magnetic field con¬ 
centrations. The authors interpreted the observed 
oscillations as primarily due to propagating slow 
magneto-acoustic waves. The energy flux carried 
by these waves was estimated to be ~ 2.5x10^ 
erg cm“^ s“^ using WKB approximation, which 
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is lower than that required for a coronal hole with 
high speed wind 8x10^ erg cni“^ 

Tian and Xia [2008] used observations from 
Transition Region and Coronal Explorer (TRACE) 
to study intensity oscillations in a network region 
at the boundary of an equatorial coronal hole. Im¬ 
ages from two different passbands 1600 A and 171 
A representative of chromosphere and low corona, 
respectively, were used in the analysis. Power maps 
constructed in three frequency regimes, high (1.2 
- 2.0 mHz), intermediate (2.6 - 4.0 mHz), and low 
(5.0 - 8.3 mHz), indicate a suppression of high fre¬ 
quency power over the magnetic network in both 
the passbands suggesting the extension of magnetic 
shadows to lower corona. Detailed investigation on 
individual features showed PDs of 5-10 min period¬ 
icity, propagating at speeds of the order of the coro¬ 
nal sound speed. Interpreting the observed PDs 
as slow magneto-acoustic waves, the authors cal¬ 
culated the energy ffux contained in these waves 
corresponding to both the passbands. For 171 A 
the value is about 40 erg cm“^ s“^, which is far 
below the energy requirement for quiet corona and 
that for 1600 A, is about 1.368x10^ erg cm“^ s“^, 
which is of the order of the required energy for 
chromosphere. 

4. DISSIPATION/DAMPING IN WAVES 

One of the most desirable characteristics of 
MHD wave observations in the solar atmosphere, 
is their dissipation. The energy carried by these 
waves has to be deposited at appropriate heights 
to facilitate coronal heating and solar wind accel¬ 
eration. It turned out that the compressive waves 
are easy to dissipate with the conventional phys¬ 
ical mechanisms whereas the incompressive waves 
require some special conditions. 

Propagating slow magneto-acoustic waves were 
known to disappear after travelling some distance 
along the guiding structure. Their amplitude de¬ 
cays rapidly as they propagate outwards. Re¬ 
sults from theoretical modelling suggest that ther¬ 
mal conduction, compressive viscosity, optically 
thin radiation, area divergence, and gravitational 
stratification affect the wave amplitude [De Moor- 
tel^ 2009, and the references therein]. Recent ob¬ 
servations by Krishna Prasad et al [2012] reveal 
that the damping in these waves is frequency- 
dependent. The authors constructed power maps 
in different periodicity ranges and found that the 
power in long-period range is significant up to 
longer distances compared to that in the short- 
period range. This dependence, evidently, comes 
from the damping mechanism and hence useful to 
identify it. To further explore this, Krishna Prasad 
et al. [2014] studied the quantitative dependence 


of damping length on oscillation period in several 
on-disk and off-limb structures. Figure 5 displays 
their results for sunspot loops and on-disk plume 
in the top panels and that for plume interplume re¬ 
gions above the limb in the bottom panel. The left 
and right panels are for two different temperature 
channels of AIA (171 A and 193 A). Clearly, the 
dependence of damping length in polar regions is 
different from that in the on-disk regions. Besides, 
neither of them seem to agree with the expected 
theoretical dependence for any of the dominant 
physical mechanisms thought to be responsible for 
damping in slow waves. Table 1 lists the expected 
theoretical dependence of damping length on os¬ 
cillation period for several physical mechanisms. 
These relations were derived following a ID linear 
MHD model [De Moortel and Hood^ 2003, 2004] 
that is applicable under the assumptions that the 
amplitude of the oscillations is small and plasma-/? 
is unity. This discrepancy between the theory and 
observations implies the current ID linear MHD 
models are inadequate in describing the damping 
in slow waves. It is possible that the waves in 
polar regions undergo non-linear steepening that 
causes enhanced viscous dissipation [Ofman et al.^ 
2000; Afanasyev and Nakariakov^ 2015] which is 
not accounted for in this model. However, the ob¬ 
servations did not indicate any signatures of non¬ 
linearity or shock formation. So it remains to be 
seen if phase mixing, mode coupling or some new 
mechanism plays a key role in slow wave damping. 
It may also be noted that more data is required 
particularly from the on-disk loops and plume¬ 
like structures to establish the exact dependence 
of damping length on the oscillation period. 

Gupta [2014] also studied the dissipation of slow 
MHD waves in a polar coronal hole above the 
limb. The author found two different dissipation 
regimes, one close to the limb (< 10 Mm) where 
the waves are heavily damped and the other at 
greater heights (10 - 70 Mm) where the waves 
damp slowly. Frequency-dependence was also in¬ 
vestigated in both these regimes by taking three 
period ranges 4-6 min, 6-15min, and 16-45 min. 
The dissipation at greater heights was found to be 
proportional to the period while that close to the 
limb do not show any preferred relation. It was 
suggested that the dissipation at greater heights 
could be due to thermal conduction but the cause 
for heavy damping close to the limb is unclear. 
A direct comparison of this study with Krishna 
Prasad et al. [2014] would not be possible since 
the results here correspond to three pre-defined fre¬ 
quency ranges whereas those from Krishna Prasad 
et al. [2014] correspond to several individual fre¬ 
quency peaks identified from the Fourier power. 

Dissipation of Afven waves by pure viscosity, re¬ 
sistivity, and thermal conductivity is not signifi- 
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Figure 5. Dependence of observed damping length on oscillation period for slow magneto-acoustic 
waves. Panels on the top display results for structures on-disk and those at the bottom display results 
for polar regions above the limb. Left and right panels are for 171 A and 193 A channels of AIA, 
respectively. Adapted from Krishna Prasad et al. [2014]. 


TABLE 1. Dependence of damping length on oscillation period (P ) for slow waves damped 
by different physical mechanisms. Prom Krishna Prasad et al. [2014] 


Physical mechanism 

Amplitude growth of 

Period dependence of 


density perturbations 

damping length {L^ oc) 

Thermal conduction 
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e ‘2dp/‘^Cs 

pO 


2 z 


Compressive viscosity 

e ^ 

p2 

Optically thin radiation 

_ ^ 

e *^5 

pO 

Gravitational stratification 


pO 

Magnetic field divergence 

Rq/t 

pO 


cant in a homogenous laminar equilibrium plasma 
at coronal temperature and density. However, the 
lower corona is not homogeneous, nor at equilib¬ 
rium, and there is ample evidence of persistent fluc¬ 
tuations and turbulence in the coronal plasma that 
may affect the dissipation of the waves. In fact, if 
the plasma is inhomogeneous, processes like phase 
mixing, resonant absorption, nonlinear dissipation. 


and mode coupling, can enhance the dissipation of 
Alfven waves. For example, Alfven waves with op¬ 
posite phases on neighboring field lines enhance the 
viscous and resistive dissipation rates [Heyvaerts 
and Priest^ 1983]. In coronal holes, perhaps the 
most relevant disssipation mechanism for Alfven 
waves is via turbulent cascades. Low-frequency 
Alfven waves get reflected due to the gradients 
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in magnetic field and density, and interact with 
the outward propagating waves causing a turbu¬ 
lent cascade that transfers energy to smaller scales. 
Upon reaching sufficiently small scales, the energy 
is dissipated to the ambient medium through wave- 
particle interactions. Several theoretical models 
based on this mechanism were proposed to ex¬ 
plain the acceleration of fast solar wind in coronal 
holes [Coleman^ 1968; Ofman^ 2004; Chandran and 
Hollweg^ 2009; Verdini et a/., 2010; Pucci et a/., 
2014]. Observationally the damping in Alfven 
waves has been detected indirectly from the vari¬ 
ations in spectral line width [Bemporad and Abbo^ 
2012; Hahn et a/., 2012]. While most of the pre¬ 
vious studies suffered from scatter light and other 
instrumental problems these authors had overcome 
some of those difficulties and seem to show a clear 
decrease in line width with height, an indication of 
Alfven wave damping. 

5. SUMMARY AND CONCLUSIONS 

MHD waves are integral part of the coronal hole 
dynamics. The uninterrupted view of full disk of 
the Sun by SDO delivering ultra high-resolution 
images and high-resolution spectroscopic observa¬ 
tions by Hinode led to a remarkable progress in our 
understanding on these waves in the past few years. 
High speed quasi-periodic upflows were observed 
in open structures and found to cause an ambigu¬ 
ity with the outward propagating slow magneto¬ 
acoustic waves. In polar plume/interplume regions 
this can be partly avoided by using the coherence 
nature of waves [Krishna Prasad et al.^ 2011] and 
it was found that the flow-like behaviour is domi¬ 
nant only near the base of the guiding structure 
[Nishizuka and Hara^ 2011; Gupta et al.^ 2012]. 
Nevertheless, it appears that these periodic up¬ 
flows at the base and slow waves are intercon¬ 
nected. Recent 3-D MHD simulations by Of man 
et al [2012] and Wang et al [2013], for a coronal 
loop, reveal that injecting tiny upflow pulses at the 
bottom inevitably excites slow magneto-acoustic 
waves that propagate along the loop, confirming 
this. The authors also indicate that small-scale 
impulsive heating events such as nanoflares at the 
loop base can produce both the upflows and waves. 
Perhaps, this could also explain the origin of long 
period (few tens of minutes) slow waves in the solar 
atmosphere. 

Besides providing the additional momentum for 
the acceleration of solar wind, the dissipation of 
slow waves, if understood well, can provide us 
with a wealth of information about the coronal 
conditions (for e.g., transport coefficients) through 
MHD seismology. The abundant observations of 
these waves in active region loops, coupled with 


the recent advances in theoretical modelling gave 
us significant understanding on their dissipation. 
Still, there appears to be a discrepancy between 
the theoretical and observed dependence of damp¬ 
ing on the oscillation period of these waves [Kr¬ 
ishna Prasad et al.^ 2014]. The observed depen¬ 
dence in plume/interplume regions was also found 
to be considerably different from that in the on- 
disk loops despite the similarity in the waves’ prop¬ 
erties. This latter disagreement could possibly 
arise from the differences in physical conditions of 
these two regions. 

The incompressive Alfven waves were believed 
to carry sufficient flux to meet the coronal hole 
energy budget, but been elusive so far. There 
were some indirect evidences from the variation 
of line broadening with height above the limb 
but most of them are crippled by the inability to 
separate thermal/non-thermal components of line 
width and contamination from scattered light had 
been another hindrance. Recent observations by 
Bemporad and Abbo [2012] and Hahn et al [2012] 
seem to overcome these difficulties. The authors 
found the width of several spectral lines initially 
increase up to certain height above the limb and 
then steadily decrease even after carefully account¬ 
ing for the scattered light. Since the ion tempera¬ 
tures do not decrease with height, this had been a 
compelling evidence for the existence and damp¬ 
ing of Alfven waves in the polar regions. The 
dissipated energy flux estimated by these authors 
also match with the requirement for coronal holes. 
However, these observations could not provide any 
information on the exact dissipation mechanism. 
The source of these oscillations also remains un¬ 
clear. We hope the recently launched Interface 
Region Imaging Spectrograph (IRIS) and the up¬ 
coming solar missions will shed more light on this 
subject. 
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